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Abstract 

Critical  conditions  for  shear  localization  in  thermoviscoplastic 
materials  are  obtained  in  closed  form  for  idealized  models  of  simple  shearing 
deformations.  The  idealizations,  which  include  the  neglect  of  heat 
conduction,  inertia,  and  elasticity,  are  viewed  as  quite  acceptable  for  many 
applications  in  which  shear  bands  occur.  Explicit  results  obtained  for  the 
idealized,  but  fully  nonlinear  problem  show  the  roles  of  strain  rate 
sensitivity,  strain  hardening,  and  initial  imperfection  on  the  localization 
behavior.  Numerical  solutions  for  two  steels  are  shown  to  exhibit  the 
principal  features  reported  for  torsional  Kolsky  bar  experiments  on  these 
steels.  Mathematically  exact  critical  conditions  obtained  for  the  fully 
nonlinear  problem  are  compared  with  critical  conditions  obtained  by  means  of 
linear  perturbation  analysis.  Use  of  relative  changes  instead  of  absolute 
changes  in  the  linear  perturbation  analysis  gives  better  agreement  with  the 
predictions  of  the  fully  nonlinear  analysis. 


1. 


Introduction 


Shear  instabilities  in  the  form  of  shear  bands  are  commonly  observed  in 
metals  and  polymers  subjected  to  large  deformations.  The  formation  of  a  shear 
band  is  often  an  immediate  precursor  to  rupture  of  the  material.  Even  when 
rupture  does  not  occur,  the  development  of  shear  bands  generally  reduces  the 
performance  of  the  material.  Thus,  improved  understanding  of  shear  band 
formation  is  critical  to  the  development  of  improved  materials  and  components 
made  from  these  materials. 

Shear  bands  can  be  divided  into  two  types:  those  in  which  thermal 
softening  plays  a  negligible  role  in  their  formation  and  those  in  which 
thermal  softening  plays  a  primary  role.  In  the  former  case  the  shear  bands, 
sometimes  called  isothermal  shear  bands,  form  as  a  result  of  strain  softening 
due,  for  example,  to  material  damage,  to  the  development  of  soft  textures,  or 
to  phase  transformations.  In  the  latter  case  the  shear  bands,  often  called 
adiabatic  shear  bands,  form  as  the  result  of  an  autocatalytic  process:  an 
increase  in  strain  rate  in  a  weaker  zone  causes  a  local  increase  in 
temperature  which  in  turn,  for  a  thermal  softening  material,  causes  a  further 
increase  in  strain  rate. 

In  this  paper  we  consider  both  types  of  shear  bands.  We  limit  our 
attention  to  simple  shearing-  deformations.  Two  fundamental  questions 
regarding  the  critical  conditions  for  shear  band  formation  are  addressed. 

1.  For  a  given  constitutive  law,  will  shear  localization  occur  for  a 
sufficiently  large  shear? 

2.  If  so,  what  is  the  critical  shear  yc,  outside  of  the  shear  band,  for 
which  the  catastrophic  process  occurs? 

As  background  for  this  study  we  note  that  an  analysis  of  the  stability  of 
homogeneous  simple  shearing  deformations  has  been  presented  by  CLIFTON  (1978) 


for  quasi-static  deformations  and  BAI  (1982)  for  dynamic  deformations.  They 
used  a  classical,  linear  perturbation  analysis  in  which  the  coefficients  in 
the  linear  differential  equations  for  the  perturbations  were  assumed  to  vary 
sufficiently  slowly  that  these  variations  could  be  neglected  in  estimating  the 
rate  of  growth  or  decay  of  fluctuations  from  the  homogeneous  solution.  This 
procedure  determines  a  critical  strain  at  which  fluctuations  begin  to  grow; 
however,  this  initial  growth  may  or  may  not  lead  to  instability  depending  on 
the  neglected  effects  of  the  time  dependence  of  the  coefficients  and  the 
nonlinearity  of  the  complete  system  of  equations.  MOLINARI  and  CLIFTON  (1983) 
and  MOLINARI  (1984,  1985)  have  presented  some  analytical  solutions  of  the 
fully  nonlinear  problem  under  quasi-static  and  adiabatic  (no  heat  conduction) 
conditions.  With  these  solutions  available  for  measuring  the  reliability  of 
more  simple  approaches  for  determining  the  onset  of  instability,  MOLINARI 
(1985),  and  FRESSENGEAS  and  MOLINARI  (1987)  developed  a  so-called  relative 
linear  perturbation  analysis  that  accounts,  in  part,  for  the  non-steadiness  of 
the  homogeneous  solution  by  linearizing  in  the  relative  perturbation  defined 
as  the  perturbation  divided  by  the  corresponding  unperturbed  quantity.  This 
approach  has  been  shown  to  give  predictions,  as  to  whether  or  not  shear  bands 
will  form,  that  are  more  in  agreement  with  the  fully  nonlinear  theory  than  are 
predictions  based  on  classical  ■  linear  perturbation  analysis.  DAFERMOS  and 
HSIAO  (1983)  obtained  a  priori  estimates  of  the  asymptotic  behavior  of  the 
solution  of  the  nonlinear  problem  (including  inertia,  but  not  heat  conduction) 
for  the  case  of  a  Newtonian  fluid  with  temperature-dependent  viscosity. 
TZAVARAS  (1984)  extended  these  results  to  the  case  of  non-Newtonian  fluids 
with  temperature-dependent  viscosities. 

Numerical  solutions  of  the  fully  nonlinear  system  of  equations  have  been 
presented  by  several  authors:  SHAWKI,  CLIFTON  and  MAJDA  (1983),  SHAWKI 


(1986),  WRIGHT  and  BATRA  (1985),  MOLINARI  (1985).  From  these  solutions  one 
can  conclude  that  dynamical  effects  and  heat  conduction  are  relatively 
unimportant  for  steel  specimens,  with  lengths  of  5-10  mm,  subjected  to 
shearing  rates  of  10^s-1  as  in  the  torsional  Kolsky  bar  experiments  of  COSTIN, 
CRISMAN,  HAWLEY  and  DUFFY  (1979),  and  HARTLEY  (1986).  Thus,  in  this  paper  we 
neglect  dynamical  effects  and  heat  conduction  in  order  to  present  an 
analytical  approach  to  the  fully  nonlinear  problem  of  thermoviscoplastic 
localization  in  simple  shear.  Our  aim  is  to  obtain  simple  analytical  formulae 
for  determining  whether  or  not  a  shear  strain  localization  instability  will 
occur  and,  if  so,  the  critical  strain  yc  at  which  the  localization  becomes 
catastrophic.  The  boundary  conditions  will,  in  some  cases,  be  general  whereas 
in  others  they  will  be  restricted  to  a  constant  imposed  shear  stress  or  a 
constant  imposed  velocity.  Isothermal  shear  bands  are  considered  in  Section  2 
and  adiabatic  shear  bands  are  considered  in  Sections  3  and  4. 


2. 


Isothermal  problem 

We  consider  a  simple  shearing  deformation  of  strain  hardening  material 
with  strain-rate  sensitivity.  For  illustration,  we  consider  the  following 
constitutive  law: 

r  =  f(7)  7m  (»  >  o)  (1) 

where  t  is  the  shear  stress,  y  is  the  shear  strain  and  y  is  the  shear  rate. 
The  function  f(7)  takes  account  of  the  strain  hardening.  This  function  is  not 
necessarily  monotonically  increasing  in  order  to  account  for  possible  strain 
softening. 

Suppose  that,  for  a  constant  applied  strain  rate  y,  the  shear  stress  r 
passes  through  a  maximum.  Will  strain  localization  occur?  By  localization  we 
mean  that  in  some  narrow  region,  the  strain  becomes  much  larger  than 
elsewhere.  More  precisely,  we  can  define  two  types  of  localization. 

Lp  localization 

If  the  heterogeneity  of  the  solution  is  growing  so  that  at  time  t,  a 
region  R  (necessarily  narrow)  exists  where  y/y  >  P  with  y  representing 
the  average  deformation  and  P  being  a  large  number,  then  Lp  localization 
of  the  deformation  is  said  to  occur. 

L„  localization 

If  for  every  point  A  different  from  B,  the  ratio  7B/7A  tends  to 
infinity  with  increasing  time,  then  L„  localization  of  the  deformation  at 
the  point  B  is  said  to  occur. 

The  analysis  of  localization  in  this  section  is  performed  in  two 
different  ways.  First  we  derive  an  analytical  solution  of  the  fully  nonlinear 
problem.  Then  an  absolute  and  a  relative  linear  perturbation  analysis 


V 


-6- 


are  performed  and  the  corresponding  predictions  are  compared  with  the  exact 
solution. 


2.1  The  nonlinear  theory 

We  consider  a  slab  with  a  geometrical  defect.  The  width  i(y)  is 
nonuniform  as  shown  in  Figure  1.  Using  the  same  approach  as  HUTCHINSON  and 
NEALE  (1977)  for  the  uniaxial  tension  of  a  bar,  we  get  from  the  equilibrium 
equation  written  at  two  different  points  A  and  B: 


*AtA  =  ak  f(7A)(7A)In  =  *BtB  =  *B  f(7B)(7B)J 


Taking  the  power  1/m  of  each  term,  we  get  after  integration: 


*A1/m  (fU)1/m  =  V/m  (f(U)1/m  <4 


where  and  yB  are  the  initial  strains  at  points  A  and  B.  If  (fU))*/m  is 
integrable  at  infinity,  then  the  values  of  the  integrals  are  finite.  While 
maintaining  the  equality  (3),  let  yB  and  he  increased  until  the  strain 
becomes  infinite  at  one  of  the  two  points,  say  B.  Then  there  exists  for  each 
a  finite  strain  for  which  Eqn.  (3)  is  satisfied.  Hence,  we  have  La, 

localization  of  strain  if  and  only  if  the  function  (fU))1//m  is  integrable  at 
infinity. 

Let  us  consider  two  examples: 

(a)  f  ( c  )-»0  as 

Assume  that  f(^)  has  power  law  behavior  at  infinity  of  the  form 


f U)  ~  arp  as 
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where  a  and  p  are  positive  constants.  Then,  from  the  integrability  condition 
(3),  the  deformation  exhibits  L.  localization  if  and  only  if 

-p  +  m  <  0  .  (5) 

This  condition  illustrates  the  stabilizing  effect  of  the  strain-rate 
sensitivity  for  m  >  0.  Even  if  the  material  is  strain  softening  (p  >  0), 
localization  will  occur  only  if  m  is  sufficiently  small  (m  <  p)  (Figure  2c). 

If  localization  does  occur,  the  critical  localization  strain  can  be 
easily  calculated.  Let  us  consider  the  following  function  h  defined  for  each 
point  M  of  the  slab: 


h(M)  =  %/m  (f(U)1/m<H.  (6) 

«yO 

Localization  of  strain  will  occur  at  the  point  B  where  the  function  h  is  a 
minimum.  The  critical  localization  strain  at  a  point  A  is  given  by  the 
implicit  equation. 


<f (^ ) )1/m  <4. 


(7) 


In  the  particular  case  where  *he-  initial  strain  y°  is  uniform  the  localization 
will  occur  at  the  points  B  where  the  width  is  minimum. 

For  m  -  p  >  0,  L*  localization  does  not  occur.  Indeed,  it  is  readily 
shown  in  Appendix  A  that 


li™  (7B/7  A)  =  (WfB)1/(m"P)  (8) 

?A~” 

This  quantity  tends  to  infinity  as  m  -  p  -  0+.  Then,  according  to  our 
definition,  Lp  localization  occurs  for  m  -  p  small  enough  (Figure  2b). 


•'  v'v'v  .  '  v; .■ 
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In  the  second  example,  we  consider  a  function  f(;)  for  which  the  limit  at 


infinity  exists 


lim  f(^)  =  L  >  0  . 

£  -•  °o 

Then  the  result  in  appendix  A  combined  with  Equation  (3)  shows  that: 

(W*B>  Vm  =  [*  fU)1/m  d;/  fU)1/1*  d;  -  (Ll/m  7b)  /  (Ll/m  7a) 

J7° 

"  7B/7A  as  7a  -  00  • 


(9) 


(10) 


The  shape  of  the  curve  7a  -  7b  is  shown  in  Figure  3  for  yB  *  As  passes 
through  the  strain  at  which  f(yB)  is  a  maximum  the  ratio  7B/7^  increases 
strongly.  It  is  then  possible  that  7B/7A  takes  on  large  values,  say  *  P. 
Thus,  Lp  localization  may  occur  just  after  f(7B)  passes  through  a  maximum.  As 
the  deformation  continues,  the  ratio  7b/7a  decreases  and  tends  to  7b/7a|°°  ” 
(iA/iB)ly/,m*  (J!A/J!B)1^m  ■'■s  not  ^ar<3e  enough,  Lp  localization  does  not  occur 
as  7  -  co.  Therefore,  Lp  localization  may  occur  after  the  maximum  of  f(7)  is 
reached,  and  disappear  for  larger  deformations. 

Strain  rate  sensitivity  has  a  strong  stabilizing  effect  on  the  asymptotic 
behavior  of  the  solution.  For  example,  consider  a  1%  geometrical  defect: 


V*A  =  °-99 


Values  of  (7B/7a)«  =  (W'b)1^"  are  given  in  Table  1  for  different  values  of  m 


Table  1: 

Influence  of  Strain  Rate  Sensitivity  on  the 
Asymptotic  Behavior  of  Plastic  Flow 
(geometrical  defect:  iB/iA  =  °*") 

m  0.2  0.1  0.01  0.006  0.001 

(7B/7A)~  1.052  1.106  2.732  7.464  2.316  X  104 


A  value  m  =  0.01  is  sufficient  to  prevent  pronounced  localization  as  y  -  ®. 
For  small  values  of  m  (say  m  <  0.005)  Lp  localization  occurs  as  y  -  °°. 

These  examples  illustrate  that  a  maximum  in  the  stress-strain  curve  does 
not  lead  necessarily  to  the  localization  of  plastic  flow. 

2.2  Linear  Perturbation  Analysis 

It  is  interesting  to  compare  the  results  of  the  fully  nonlinear  theory  to 
the  predictions  of  a  linear  stability  analysis.  For  this  comparison  consider  a 
block  of  uniform  thickness  *(y)  =  i0  undergoing  homogeneous  simple  shearing 
deformation  70(t).  Let  67  =  7(y,t)  -  70(t)  be  the  difference  between  the 
shear  strain  7(y,t)  for  the  same  block  subjected  to  the  same  boundary 
conditions,  but  having  a  fluctuation  in  strain  and  strain  rate  beginning  at 
some  time  t0.  Using  the  constitutive  law  (1)  and  considering  the  problem  as 
quasistatic  (6  t  =  o)  we  obtain: 


67  decays  initially.  Because  the  right  side  of  Equation  (9)  is  independent  of 
the  coordinate  y,  the  growth  or  decay  of  the  strain  difference  is  such  that 
the  ratio  of  the  strains  syA  and  573  at  points  A  and  B  remains  constant.  Thus, 
Equation  (9)  provides  no  information  on  strain  localization. 

From  this  linear  analysis  we  see  that  whether  small  perturbations  of  a 
homogeneous  shearing  deformation  are  expected  to  grow  or  decay  initially 
depends  only  on  the  sign  of  f'(70),  the  slope  of  the  stress-strain  curve  at  a 
constant  strain  y0.  However,  from  the  nonlinear  analysis,  we  know  that 
whether  or  not  localization  will  occur  is  not  governed  by  the  value  of  f'(70), 
but  by  the  strain  rate  sensitivity  parameter  m  and  the  behavior  of  f(y)  as 
7-“>.  The  condition  f'(7)  <  0  is  a  necessary  and  sufficient  condition  for  the 
initial  growth  of  perturbations  of  a  homogeneous  deformation  whereas  f'(7)  <  0 
for  7  greater  than  some  critical  value  yc  is  only  a  necessary  condition  for 
localization  to  occur.  This  tendency  for  the  linear  perturbation  analysis 
based  on  (9)  to  predict  the  growth  of  small  perturbations  under  relatively 
weak  restrictions  on  the  constitutive  equations  can  be  partially  offset  by 
considering  the  relative  perturbation 
67 

Ay  =  —  .  (10) 

?o 

The  relative  perturbation  A7  tends  to  grow  more  slowly  than  the  absolute 
perturbation  67  and  may  even  decay  as  the  perturbation  grows. 

From  (9)  and  (10)  we  obtain  after  logarithmic  differentiation 


Integration  of  (11)  gives 


f(70)"1/m 

A7  =  K  - 

7o 


(12) 


If  f(y0)  has  the  behavior  (4)  for  large  values  of  yQI  then  as  70-®  the 
relative  perturbation  A7  becomes  unbounded  for  -p  +  m  <  0  and  approaches  zero 
for  -p  +  m  >  0.  These  conditions  are,  respectively,  the  same  as  the  critical 
conditions  (See  example  (a)  of  the  previous  section)  for  Lp-localization  to 
occur  or  not.  This  parallellism  between  predictions  of  the  linear  relative 
perturbation  analysis  and  the  exact  results  for  the  nonlinear  theory  suggests 
that  linear  relative  perturbation  analysis  may  be  more  widely  useful  in 
predicting  the  stability  of  deformations  than  is  the  commonly  used  linear 
perturbation  analysis  represented  by  Equation  (9).  However,  we  emphasize  that 
the  localization  analysis  in  the  nonlinear  theory  and  the  linear  relative 
perturbation  analysis  address  different  problems  and  there  is  no  a  priori 
reason  to  expect  that  the  critical  conditions  for  L^-localization  are,  in 
general,  the  same  as  the  critical  conditions  for  predicted  unbounded  growth  of 
a  relative  perturbation. 
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3. 


Adiabatic  Case 


We  consider  next  the  influence  of  temperature  on  localization.  As 
discussed  in  the  introduction,  the  deformation  is  assumed  to  be  adiabatic  and 
quasistatic.  We  consider  the  constitutive  equation 

r  =  t  (7,  7,  0),  (13) 

the  equation  of  equilibrium 

i(Y)  T(y,t)  =  i(h)  T(h,t) ,  (14) 

che  compatibility  equation 
3v 

y  =  —  *  (is) 

3y 

and  the  energy  equation 


30 

pC  —  =  6  t  7  .  (16) 

at 

In  these  equations  p  is  the  mass  density,  C  is  the  heat  capacity  per  unit 
mass,  0  is  the  absolute  temperature,  v  is  the  particle  velocity,  and  0  is  the 
Taylor-Quinney  coefficient  which  characterizes  the  fraction  of  plastic  work 
that  is. converted  into  heat;  usually  0  is  taken  constant  and  equal  to  0.9. 
Equations  (13)-(16)  constitute  four  equations  in  the  four  unknowns  7,0,r,v. 

In  the  following,  we  present  a  discussion  of  localization  for  different 
constitutive  laws  and  different  boundary  conditions.  We  consider  the  cases  of 
constant  velocity  boundary  conditions: 


v(0,t)  =  0 
v(h#t)  =  v0 

or  constant  stress  boundary  conditions 


jB(h)T(h,t)  =  i(0)r(0,t)  =  const. 


3.1  Materials  Without  Strain-Hardening 

An  exact  solution  of  the  fully  nonlinear  problem  has  been  presented  by 
MOLINARI  and  CLIFTON  (1983)  for  the  case  in  which  the  material  is  not  strain 
hardening  and  Eqn.  (13)  has  the  form 

t  =  4(0)  ym  .  (19) 

In  order  to  obtain  this  exact  solution  we  write  Eqns.  (14)  and  (16)  at  two 
different  points  A  and  B.  Substitutions  of  (19)  into  (14)  and  use  of  (16) 
to  eliminate  7^/73  gives 


4m+1)/mM(eA)1/m  d0A  «  iBm+1)/I\(0B)1/in  d0j 


which,  after  integration,  becomes 


(m  +  l)/m  [  A 


M(5)1/m  dt  =  *B<m  +  i)/®  M(C)1/m  cH 


where  0A  and  8§  are  the  initial  temperatures  at  points  A  and  B. 

From  (21)  it  appears  that  L*.  localization  of  temperature  occurs  at  B  if, 
and  only  if,  //(G)1/111  is  integrable  at  infinity,  i.e.  if  there  exists  some  K  > 
0  for  which 
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I 

*  *  - 

CV. 

I 

to 

y 

I 


M(5)1/m  <  +  •  . 


(22) 


K 


Localization  will  occur  at  the  point  B  where  the  following  function,  defined 
for  each  point  M  of  the  slab, 


M  -  iM(m  +  l>/» 


M(U1/m  <4 


(23) 


a° 

eM 


is  a  minimum.  At  localization  the  temperature  0^  at  any  point  A  is  given 
by 


f  e 


i  (m  +  l)/m 
A 


n(  5 )l/m  <4  =  i(m  +  1)/m 

B 


mU)1/®  <4  . 


(24) 


0 


B 


It  is  easy  to  show  that  L,*  temperature  localization  tends  to  result  in 
strain  localization.  Indeed,  from  the  equilibrium  condition  (15)  and  the 
constitutive  law  (19),  we  have 


*A  ^(0A)  7l  =  <B  ^<0b)7b 

Then,  assuming  temperature  localization,  we  get 

,1/m 


,m 


lim 

'  ?B  ' 

_  lim 

*h  *40a) 

eA  -  eA 

.?A  . 

eA  ■*  eA 

.  *b  *40b)  . 

=  +  oo 


since  lim  4(83)  =  0  from  t^e  integrability  condition  (22).  This  L« 


0A  -  0A 


localization  of  the  strain  rate  essentially  ensures  L„  localization  of  the 
strain  although  various  pathological  cases  must  be  excluded  in  a  rigorous 
analysis.  We  henceforth  consider  constitutive  equations  and  loading  conditions 


for  which  such  pathological  cases  are  excluded. 

Localization  results  obtained  from  the  integrability  condition  (22)  are 
summarized  in  Table  2  for  several  constitutive  laws. 

Table  2:  Localization  results  for  visco-plastic,  thermal  softening 
materials  without  strain  hardening 
(m  >  0,  40  >  0,  a  >  0) 


Constitutive  law 

L*  localization 

L1 

r  =  n0  ev  7m 

v  +  m  <  0 

L2 

t  =  u0  exp(-«9)  r® 

a  >  0 

l3 

t  =  sup (a  +  b0,O)  7® 

b  <  0 

L4 

t  =  n0  exp( -a/0 )  7® 

never  exhibits 

L„  localization 
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3.1.1 


Calculation  of  the  critical  strain 


In  this  section  we  obtain  explicit  results  for  the  critical  strain  at 
localization.  We  illustrate  the  approach  by  considering  the  constitutive  law 
L2.  Substituting  this  law  in  Eqn.  (21),  we  obtain  the  following  expression  for 
the  temperature  at  a  point  B  as  a  function  of  the  temperature  at  a  point  A 


0B  =  -  —  log  I  — 


ja'|(1  +  m)/m 


exp(-a0A/m)  +  Cj 


where 


C1  =  exp(-a9g/m)  -  — 


i A  ■)  ( 1  +  n)/m 


exp(-a0A/m) 


From  (25),  a  necessary  condition  for  localization  to  occur  is  <  0 
since  for  Cj  ?  0  the  logarithm  cannot  tend  to  infinity.  Let  us  identify  the 
point  B  as  the  point  where  the  quantity  jB(ln+1)/111  exp(-a0g/m)  is  a  minimum. 
Then,  if  the  initial  temperature  (or  the  width  i)  is  non  uniform,  ^  is 
strictly  negative. 


For  a  thermal  softening  material,  i,e.  a  >  0,  the  quantity  expf-crO^m) 
decreases  to  zero  as  0A  -  ».  For  some  critical  temperature  0A,  the 
temperature  0B  will  become  infinite.  From  (26),  this  value  is 


/■»  1U  LJ 

»A  -  »A  ■  i  1»9  1  - 


^1  +  m 

'B  ]  m 


exp(-a(0g  -  0g)/m) 


If  the  material  is  thermal  hardening,  i.e.  a  <  0,  then  the  term  exp(-a0A/m) 


grows  and,  from  (25),  it  is  obvious  that  L*  localization  is  impossible. 


Indeed,  for  sufficiently  large  0A,  the  difference 

rr  ih  i  -I 

- l0?r  [( -  J  ] 

becomes  small  compared  to  the  absolute  temperature,  say  0A. 

.  To  calculate  the  critical  strain  from  the  critical  temperature  we 

consider  first  the  case  in  which  a  constant  stress  r  is  applied  at  the 
boundary.  Then,  from  the  energy  equation  (16),  0  can  be  calculated  as  a 
function  of  the  strain  y 
Bt 

9  =  —  7  +  e°.  (28) 

pc 

The  critical  strain,  at  A  is  obtained  by  substituting  the  critical 

temperature  0^  given  by  (27)  into  (28)  to  obtain 

1  +  m 

pCm  flg 

7a  =  •  -  l0<?  1  “  —  exp(-a(0g  -  0g  )/m)  .  (29) 

a0T  L  l  JA 

The  stabilizing  effects  of  increased  strain  rate  sensitivity  (i.e.  larger  m) 
and  decreased  thermal  softening  (i.e.  smaller  a)  are  evident  in  this 

expression.  The  relative  importance  of  geometrical  and  temperature  defects  is 
evident  from  (29)  which  shows  that  the  critical  strain  yA  is  the  same  when 
(ib/«a)  an<i  ( ®b  ”  ®a)  to  the  locus  of  values  for  which  the  argument  of 

the  logarithm  in  (29)  is  a  constant.  For  small  values  of 
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This  approximation  is  discussed  in  appendix  B  and  a  numerical  evaluation  will 
be  presented  later.  Substitution  of  the  approximate  stress  (34)  into  the 
energy  equation  (16)  leads  to 


i  30 

;  PC  —  =  0  m0  exp(-«9)  Jo  7  • 

)  .  at 


(35) 


This  equation  can  be  integrated  by  separation  of  the  variables  0  and  y  to  give 


0(7)  =  0°  +  -  log 


,m 


1  + 


«  6  7S 

pC 


exp(-a0°)  7 


(36) 


where  0°  is  the  initial  temperature. 

With  the  relationship  between  0  and  y  given  by  (36),  the  critical  strain, 
7^,  at  A  can  be  obtained  by  integration  of  the  equilibrium  equation  (14).  Such 
integration  gives 


-  1/m 


1/m 

exp(-«0A(7)/m)  dy  =  iB' 


/B 


exp(-a0B(7)/m)  dy 


(37) 


0  o 

with  9a(7 )  and  0B(7)  given  by  (36).  At  localization,  yB  becomes  infinite  and 
the  critical  strain  yA  at  point  A  becomes,  for  m  «  1, 


~yC  _ 

7 A 


pc 


a  Bt9 


A  J 


1/m 


-m/(l-m) 


exp(-a(l-m)(0g  -  e£)/m) 


-  1 


(38) 


where 

tA  =  H  ?o  exp(-«0A) 


is  the  shear  stress  at  A  in  an  isothermal  deformation  at  the  same  strain  rate. 
If  I  and  0°  are  both  uniform,  then  Eqn.  (38)  implies  that  7^  is  infinite  and 
localization  does  not  occur.  The  critical  strain  decreases  as  (iB/*A) 


decreases  and  8g  -eg  increases.  The  energy  measure,  t£  yg,  of  the  critical 
strain  increases  with  increasing  strain  rate  sensitivity  (i.e.  increasing  m) 
and  decreasing  thermal  softening  (i.e.  decreasing  a).  Again,  the  nominal 
critical  strain  is  obtained  by  the  substitution  of  (38)  into  (32).  Comparison 
of  (38)  and  (29)  indicates  that,  for  m  <<  1,  the  relationship  (31)  between 
equivalent  temperature  defects  A0  and  geometric  defects  Ai  holds  for  velocity 
boundary  conditions  as  well  as  for  constant  stress  boundary  conditions. 

Further  understanding  of  the  dependence  of  the  critical  strain  yg  on  the 
defect  can  be  obtained  by  introducing  the  defect  parameter 

4  g 

X  =  1  -  _  exp  (-«(0§  -  eg)).  (39) 

L  'A  J 

In  terms  of  x  the  critical  strain,  for  m  «  1,  is  given  by 

pc  f 

7ft  -  - —  U  -  (1  -  Xl1/®]-®  -  1  |  .  (38)' 

«<»A  l 


For  sufficiently  small  defects  (i.e.  0  <  x  «  m  «  1)  this  expression  can  be 
approximated  by 


-m 


or,  alternatively, 


yg  »  -  [-m  log  x  +  const.]  (40b) 

aerg 


where  the  constant  is  chosen  such  that  (38)  and  (40b)  give  the  same  value  for 


7A  for  one  value  of  x  in  the  interval  of  interest.  From  (40b)  it  is  evident 
that  the  critical  strain  depends  weakly  on  the  defect  parameter  for 
sufficiently  small  defects.  For  values  of  m  that  are  characteristic  of  steels 
(e.g.  0.01  «  m  <;  0.02)  the  approximate  relation  (40b)  provides  a  satisfactory 
representation  of  the  dependence  of  7A  on  x  over  three  decades  of  variation  of 
X  for  x  <  0.2  m. 

Identical  calculations  can  be  performed  for  a  power  law  dependence  of  the 
flow  stress  on  the  temperature.  Analogous  results  for  the  constitutive  law 
Lj  of  Table  2  are 


1  +  Tft  V  +  Tft  m 

c  -  r\  ^  _ 


c  pC0A 

'  *B  ' 

m 

eB 

m 

v  +  m 

= - 

1  - 

— 

— 

-  1 

6rA 

.  *A  . 

l6a  J 

. 

- 

(41a) 


for  constant  stress  boundary  conditions,  and 


PC9A 

(l-v)Br^ 


v+m(l-v) 


m 


(l-v)m 
v+m( 1-v) 

-  1 


(41b) 


for  constant  velocity  boundary  conditions;  Eqn.  (41b)  holds  only  for  m  <<  1. 


3 . 2  Materials  with  Strain  Hardening 

Strain  hardening  cannot  be  ignored  for  most  materials.  In  this  section 
we  derive  analytical  localization  criteria  for  constitutive  laws  of  the  form 

t  =  n(Q)  (7  +  7°)n  7m  (42) 


where  70  is  the  initial  strain.  The  approach  is  similar  to  that  used  in 
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3.2.1  Stress  boundary  condition 

Elimination  of  y  between  the  constitutive  law  (42)  and  the  energy 
equation  (16)  leads  to 


m  +  1 

d8  Br  ^  -1/m 

—  =  —  m(9) 

dt  pC 


(7  +  7°) 


We  write  (43)  at  two  different  points  A  and  B,  take  account  of  the  equilibrium 
equation  (14),  use  (28)  to  eliminate  y,  and  integrate  to  obtain 


m  +  1  0; 


MU)  U-0X)  —  +  7a 


m  +  1 


M(t)  (5“0§) 


Equation  (44)  is  a  generalization  of  (21)  to  strain  hardening  materials. 
Analysis  of  (44)  analogous  to  that  of  (21),  shows  that  L*,  localization  occurs 
if  and  only  if  the  function 


1/m  pC  n/m 

0  -  n(B)  (6  -  0°)  —  +  7° 

8t 


is  integrable  at  infinity.  For  the  constitutive  law 


t  =  Mi  6V  (7  +  7°)n  7m 
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Lco  localization  occurs  if  and  only  if 

v  +  n  +  m  <  0  .  (46) 

The  inequality  (46)  provides  a  good  illustration  of  the  competition  between 
the  stabilizing  effects  of  strain  hardening  (n  >  0)  and  positive  strain  rate 
sensitivity  (m  >  0),  and  the  destabilizing  effects  of  thermal  softening 
(v  <  0).  The  localization  criteria  (46),  obtained  by  MOLINARI  and  CLIFTON 
(  1983  ),  has  also  been  obtained  by  FRESSENGEAS  and  MOLINARI  (1986)  as  the 
criterion  for  the  initial  growth  of  a  fluctuation  based  on  a  linear  relative 
perturbation  analysis.  The  inequality  (46)  differs  from  the  condition 

v  +  n  <  0  (47) 

that  must  be  satisfied  for  the  initial  growth  of  a  fluctuation  according  to 
absolute  linear  perturbation  analysis.  The  difference  between  the  conditions 
(46)  and  (47)  illustrates  the  tendency  for  absolute  linear  perturbation 
analysis  to  predict  growth  of  fluctuations  under  some  conditions  for  which  the 
full  nonlinear  analysis  predicts  that  localization  will  not  occur. 

For  the  constitutive  law 

T  =  //0  e_ot0(7  +  y°)n  7m  (48) 

a  similar  analysis  shows  that  L*  localization  occurs  if  and  only  if  a  >  0 
(thermal  softening).  For  this  constitutive  equation  the  critical  strain  7^  at 
A  when  the  strain  at  point  B  becomes  infinite  is  obtained  by  the  substitution 
of  (48)  into  (44)  to  obtain,  after  a  change  of  variable, 
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>_v  vn/r-  dv  =  Kr 


dv 


where 


Jh  =  —  (7A  +  7h)  >  ?A 

m  pC 


vo 
-  ?A 


1  +  n  +  m 


KA  =  lA 


with  A  replaced  by  B  for  and  Kg.  Localization  will  occur  at  the  point  B, 
where  the  quantity  on  the  right  side  of  (49)  is  a  minimum. 

Equation  (49)  has  the  same  form  as  that  obtained  by  HUTCHINSON  and  NEALE 
(1977)  in  the  study  of  the  rupture  of  a  viscoplastic  bar  in  tension  although 
the  physical  effects  being  modeled  are  different  -  their  analysis  included 
necking,  but  did  not  include  the  thermal  softening  which  is  included  here. 


3.2.2  Velocity  boundary  conditions 

As  in  section  3.1  we  consider  constant  velocity  boundary  conditions  and 
assume  that  the  strain  rate  sensitivity  of  the  material  is  small 
(i.e.  m  «  1).  To  calculate  the  temperature  from  the  energy  equation  (16)  we 
replace  y  by  y0  =  V/h  in  the  constitutive  equation  (42)  and  integrate  to  obtain 


1  B  7q 

-  =  -  (7  +  7°)n+l  -  (7°)n+l  > 

nil)  PC  (n  +  1) 


Substitution  into  (50)  of  functions  * i(9)  that  model  the  temperature 
dependence  of  the  flow  stress  gives  the  required  relationship  between  the 
temperature  0  and  the  shear  strain  7.  For  p(0)  =  ^10v  we  obtain 


0(7)  =  0° 


1  +  (1  “  v)  C  7° _  [(7  +  7°)n+1  -  (7°)n+1] 

pC  (n  +  l)(0°)1-v 


1-v 


(51) 


and  for  p(0)  =  p0  e"a9  we  obtain 


0(7)  =  0°  +  -  log 

a 


m  a-a 


1  + 


“No  ?o  e 


0° 


pC  (n  +  1) 


(7  +  yOjn+1  _(70)  n+1 


(52] 


These  equations  provide  an  approximate  relationship  between  the  temperature 
and  the  strain  at  each  position  as  long  as  the  exponent  m  is  sufficiently 
small  for  the  dependence  of  the  shear  stress  on  strain  rate  to  be  represented 
by  7q,  where  y0  is  the  nominal  strain  rate,  instead  of  by  7m,  where  7  is  the 
local  strain  rate. 

In  order  to  investigate  the  critical  conditions  for  localization,  we 
substitute  the  functions  0(7)  obtained  from  (51)  or  (52)  into  the  equation 


1/m 

*A 


’  ?A 

M(eAU))1/m  U  +  7&}n/m 


(53) 


1/m 

*B 


>  ?B 

4(6BU))1/m 


U  +  7§)n/m 
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As  before,  L*,  localization  occurs  if,  and  only  if,  the  integral  on  the  right 
side  of  (53)  remains  bounded  as  73-“.  For  8(7)  given  by  Eqn.  (51),  the 


condition  for  L*  localization  becomes 
v  +  n  +  m  (1-v)  <  0 


for  0  <  m  «  1,  and  v  <  1.  This  condition  is  slightly  more  restrictive  than 
the  condition  (46)  obtained  for  stress  boundary  conditions.  That  is,  the 
tendency  for  localization  is  slightly  stronger  for  stress  boundary  conditions 
than  for  velocity  boundary  conditions  in  that  the  localization  condition  (46) 
is  satisfied  by  all  v,m,  and  n  which  satisfy  (54);  however,  for  m  «  1,  the 
terms  involving  m  in  both  (46)  and  (54)  are  often  so  small  that,  effectively, 
the  localization  conditions  (46)  and  (54)  are  the  same.  For  0(7)  given  by 
Eqn.  (52)  the  condition  for  L*,  localization  is  satisfied  for  all  «  >  0 
provided  that  m,  n  satisfy  m  >  0,  n  >  -  1. 
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4 .  Numerical  Example 

Dynamic  torsion  experiments  for  investigating  shear  localization  have 
been  performed  by  HARTLEY  et  al.  (1986)  on  two  different  types  of  steel:  CRS 
1018  and  HRS  1020.  At  the  strain  rates  (I03s'1)  and  temperatures  (6°  =  300°K) 
of  these  experiments  the  behavior  of  these  materials  can  be  represented 
reasonably  well  by  a  constitutive  equation  of  the  form  (45).  Numerical  values 
of  the  various  parameters  in  the  model  are  given  in  Table  3  (SHAWKI  (1986)). 
The  strain  y°  is  taken  to  have  the  value  0.01  for  both  steels.  More  detailed 
fitting  of  the  plastic  response  of  these  steels  has  been  presented  by 
KLEPACZKO  (1986). 

Table  3  Thermomechanical  Properties  of  CRS  1018  and  HRS  1020  Steels 


Steel 


CRS  1018 


HRS  1020 


Parameter 


V 

-  0.38 

-  0.51 

n 

0.015 

0.12 

m 

0.019 

0.0133 

P 

7800  kg/m3 

7800  kg/m3 

C 

500  J/kg°  K 

500  J/kg°  K 

i 

3579  x  106  S.I. 

7587  x  106  S.I. 

Variations  i( y)  in  the  wall  thickness  of  the  specimens  were  not  reported 
by  HARTLEY  (1986).  Subsequently,  DUFFY  (1986)  has  sectioned  specimens  used  in 
such  experiments  to  determine  the  variation  in  wall  thickness,  both  along  the 
length  of  the  specimen  and  around  its  circumference.  For  CRS  1018  the  wall 
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thickness  is  relatively  uniform  around  the  circumference,  but  strong  variation 
-  up  to  10%  -  occur  along  the  length  of  the  specimen.  For  the  purpose  of  this 
numerical  example  we  take  the  geometrical  factor  (^b/j!a)  *n  the  Preceding 
analysis  to  be  a  parameter  that  varies  from  0.9  to  0.999.  In  order  to  relate 
the  critical  strains  to  the  nominal  strain  yc  at  localization  (see  Eqn. 
(32)),  the  variation  in  wall  thickness  iA  =  s( yA)  must  be  prescribed  over  the 
entire  length  of  the  specimen.  Based  on  the  general  appearance  of  the 
sectioned  specimens  we  take  this  variation  to  have  the  form 
i(y)  «  2ny 

-  =  1  +  -  (cos  -  -  1)  (55) 

*A  2  h 

where  £  is  a  geometrical  parameter  that  is  taken  to  vary  from  10_1  to  10~6  to 
give  the  range  of  values  of  0.9  to  0.999999  for  fB/iA. 

Boundary  conditions  for  the  dynamic  torsion  ("torsional  Kolsky  bar") 
experiment  are  effectively  those  of  imposed  constant  velocity  at  the  ends  of 
the  specimen.  Hence,  we  use  the  solution  for  velocity  boundary  conditions 
given  by  Eqns.  (51)  and  (53).  The  restriction  to  m«l  that  is  required  in 
obtaining  (51)  is  well  satisfied  by  the  values  m  =  0.019  for  the  CRS  and 
m  =  0.0133  for  the  HRS.  Evaluation  of  7c(yA)  from  (53)  and  integration  over 
the  length  of  the  specimen,  according  to  (32)  gives  the  dependence  of  the 
critical  strain  yc  on  the  geometrical  imperfection  parameter  £  that  is  shown 
in  Figs.  4  and  5.  For  small  €  the  nominal  critical  strain  varies 
approximately  as  log  €,  as  predicted  for  the  local  critical  strain  by  Eqn. 
(40b).  In  each  figure  the  insert  provides  an  expanded  scale  of  the  region  of 
primary  interest  in  the  interpretation  of  torsional  Kolsky  bar  experiments. 
For  one  value  of  €  (£  =  0.02),  the  strain  distribution  at  localization  for  CRS 
1018  is  shown  in  Fig.  6.  The  width  of  the  band  of  intense  shear  (say,  the 


region  for  which  7(y)  >  3y ( 0 ) )  is  approximately  20%  of  the  length  of  the 
specimen.  Such  relatively  wide  bands  are  observed  in  CRS  1018  (HARTLEY  et  al. 
(1985)).  The  predicted  nominal  stress-strain  curves  for  the  two  steels,  with 
the  geometrical  imperfection  parameter  t  equal  to  0.02,  are  shown  in  Figs.  7 
and  8;  a  corresponding  curve  for  e  =  0.04  is  included  in  Fig.  8.  The 
general  features  of  the  curves  include  a  slowly  rising  segment  during  which 
the  shearing  is  quite  uniform,  a  slowly  falling  segment  during  which  a  broad 
band  of  enhanced  shearing  develops,  and  a  sharply  falling  segment  during  which 
the  shearing  becomes  intensely  localized  in  a  band.  These  general  features  are 
characteristic  of  the  experimental  records  obtained  in  such  experiments 
(HARTLEY,  et  al.  (1986)).  Numerical  values  for  the  strain  at  the  peak  of  the 
stress-strain  curve  and  the  strain  at  the  beginning  of  the  sharp  decline  in 
stress  are  comparable  to  values  obtained  in  experiments.  However,  the 
predicted  rate  of  sharp  decline  is  greater  than  normally  measured.  This  rate 
of  decline  is  affected  by  the  detailed  geometry  of  the  initial  imperfection 
which  probably  was  not  modeled  adequately  by  the  generic  form  (55)  us  d  to 
model  the  imperfection.  Other  difficulties  with  comparisons  between  theory 
and  experiment  in  the  steeply  falling  part  of  the  curve  include:  (i)  the 
inadequacy  of  the  assumption  that  the  stress  obtained  using  the  nominal  strain 
rate  can  be  used  in  calculating  the  local  rate  of  energy  dissipation,  (ii)  the 
likelihood  that  the  final  localization  varies  so  strongly  around  the 
circumference  of  the  specimen  that  a  one-dimensional  analysis  is 
inappropriate;  and  (iii)  the  lack  of  constant  velocity  boundary  conditions 
when  the  stress  decreases  strongly  in  torsional  Kolsky  bar  experiments. 
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5.  Conclusion 

By  assuming  the  deformation  to  be  adiabatic  and  quasi-static,  and  by 
neglecting  elasticity  effects,  we  have  characterized,  analytically,  the 
critical  conditions  for  shear  strain  localization  in  simple  shear.  The 
assumed  conditions  are  good  approximations  for  the  specimen  sizes  and  strain 
rates  that  are  commonly  used  in  torsional  Kolsky  bar  experiments  on  shear  band 
formation  in  steels. 

We  assume  the  existence  of  initial  inhomogeneities  which  are  either 
geometrical  defects  or  non-uniform  fields  of  initial  temperature  or  strain. 
The  localization  strain  is  obtained  as  a  function  of  these  defects,  the 
material  parameters  and  the  boundary  conditions.  Two  types  of  boundary 
conditions  have  been  considered: 

-  constant  applied  stress 

-  constant  applied  velocity; 

in  the  latter  case,  the  analytical  results  are  restricted  to  materials  with 
weak  strain  rate  sensitivity. 

The  results  are  particularly  simple  for  materials  without  strain 
hardening.  In  this  case,  explicit  expressions  are  obtained  for  the  dependence 
of  the  critical  strain  on  a  defect  parameter  that  characterizes  the 
geometrical  defect  and  the  nonuniformity  of  the  initial  temperature.  For 
materials  with  weak  strain  rate  sensitivity  the  critical  strain  depends  weakly 
(essentially  logarithmically)  on  the  amplitude  of  the  imperfection  for  small 
imperfections . 

Comparison  of  predictions  of  the  theory  with  experimental  results  for  a 
cold-rolled  steel  shows  good  agreement  in  the  qualitative  features  of  the 
response.  Quantitative  comparisons  require  detailed  descriptions  of  the 
geometrical  defects  of  the  specimens  used  in  the  experiments.  Preliminary 
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comparisons  based  on  approximate  representations  of  the  geometrical 
imperfections  of  the  specimens  suggest  that  good  quantitative  agreement  may  be 
obtained  once  the  defects  are  modeled  accurately. 
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Consider  two  piece-wise  continuous  functions  g  and  h.  The  function  g(x) 
is  positive  and  nonintegrable  at  infinity, i.e.  g(;)  d£  =  +  «  .  If  the 
function  h(;}  is  equivalent  to  g(£)  at  infinity  (i.e.  lim  h(;)/g(;)  =  1),  then 

£  — *oo 

(DIEUDONNE  (1968)) 

x  x 

h(0  d;  “  g(;)  d;  ,  as  x  -  ®  .  (Al) 

a  a 


With  the  choices  h(;)  =  ff;)1/131  and  g(;)  =  a1/131  £~P/m  we  have,  from  Eqn.  (4), 
h  ~  g.  Then  Eqns.  (3)  and  (Al)  imply  that  for  m  -  p  >  0  we  have 
m  .  m 


jjl/m  al/m 


7A(m-p)/m  -,1/m  al/m  _  ^(m-pj/m 

m-p  m-p 


(A2) 


since,  when  -yA  tends  to  infinity,  must  also  tend  to  infinity.  From  (A2)  it 
follows  that 


lim  (7B/7A)  =  (ViB)1/(II'"P)  • 

7a"°° 


(A3) 


Appendix  B 


During  the  localization  process  the  local  strain  rate  7  will  differ  from 
the  nominal  strain  rate  yQ  =  V/h.  As  a  measurement  of  this  difference  we 
define  the  ratio 

,  X  =  7/70  .  (Bl) 

We  want  to  evaluate  the  acceptability  of  the  approximation  of  replacing  the 
stress 

r  -  n  e~a  9  7m  (B2 ) 

by  the  quantity 

t0  =  fi  e’«  9  7q  (B3 ) 

in  calculating  the  heat  generated  by  plastic  working.  For  m  =  0.01,  the  ratio 
t /t o  =  xm  is  bounded  by  0.955  «  t/t0  «  1.071  for  x  in  the  interval  [10-3, 
lo3].  Thus,  using  the  approximate  stress  t0  leads  to  a  maximum  error  of  7% 
for  a  variation  in  strain  rate  of  six  orders  or  magnitude. 
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